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Abstract: Recently, the ATLAS and CMS detectors have discovered a bosonic parti- 
cle which, to a reasonable degree of statistical uncertainty, fits the profile of the Standard 
Model Higgs. One obvious implication is that models which predict a significant departure 
from Standard Model phenomenology, such as large exotic (e.g., invisible) Higgs decay 
or mixing with a hidden sector scalar, are already ruled out. This observation threatens 
the viability of electroweak baryogenesis, which favors, for example, a lighter Higgs and 
a Higgs coupled to or mixed with light scalars. To assess the broad impact of these con- 
straints, we propose a scheme for classifying models of the electroweak phase transition 
and impose constraints on a class-by-class basis. We find that models, such as the MSSM, 
which rely on thermal loop effects are severely constrained by the measurement of a 125 
GeV Higgs. Models which rely on tree-level effects from a light singlet are also restricted 
by invisible decay and mixing constraints. Moreover, we find that the parametric region fa- 
vored by electroweak baryogenesis often coincides with an enhanced symmetry point with 
a distinctive phenomenological character. In particular, enhancements arising through an 
approximate continuous symmetry are phenomenologically disfavored, in contrast with 
enhancements from discrete symmetries. We also comment on the excess of diphoton 
events observed by ATLAS and CMS. We note that although Higgs portal models can 
accommodate both enhanced diphoton decay and a strongly first order electroweak phase 
transition, the former favors a negative Higgs portal coupling whereas the latter favors a 
positive one, and therefore these two constraints are at tension with one another. 
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1 Introduction 

A number of baryogenesis mechanisms are capable of explaining the observed baryon 
asymmetry of the universe, but many of these operate at a high scale, inaccessible to direct 
laboratory tests, where they evade independent confirmation. The primary motivation for 
studying electroweak baryogenesis (EWBG) [ ] is that the baryon asymmetry is generated 
by electroweak scale physics, which is tested by experiments aimed at understanding the 
nature of electroweak symmetry breaking. These include Higgs searches at LEP, the Teva- 
tron, and the LHC colliders. Thus, models of the electroweak sector may be constrained 
from two sides: by the requirement that electroweak baryogenesis successfully generates 
the baryon asymmetry and by the requirement that models remain consistent with Higgs 
search constraints. Indeed, the ATLAS and CMS collaborations recently announced the 
discovery of a particle in the mass range 125 — 126 GeV which matches the profile of the 
Higgs boson [ , ]. Even at this early stage, without a precise knowledge of the alleged 
Higgs' couplings to Standard Model (SM) fields, we have gained a partial picture of the 
origin of the electroweak symmetry breaking (EWSB). In this paper, we would like to 
understand what is the main implications of a 125 GeV SM-like Higgs for electroweak 
baryogenesis. 

Studies of the viability of electroweak baryogenesis and the impact of collider con- 
straints are usually performed on a model-by-model basis. However, many individual 
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models can accommodate a partial picture of the electroweak symmetry breaking sector. 
Thus, as the LHC begins to expose the Higgs sector, revealing only glimpses of the full 
picture, one would like to understand what classes of models may be consistent with or 
at tension with the data. To this end, we propose a scheme for classifying models of the 
electroweak sector based upon the nature of the electroweak phase transition (EWPT) and 
study the implications of the recent Higgs discovery at the LHC on a class-by-class basis. 
We find, the LHC's detection of a 125 GeV Higgs in conjunction with constraints on ex- 
otic decay and hidden sector mixing, provide strong constraints on certain EWPT model 
classes.^ 

We identify the phase transition model classes in the following way. The success of 
EWBG relies upon the electroweak phase transition being of the first order [1]. In the 
context of the phase transition calculation, this translates into the requirement that the 
thermal effective potential, K;fr(/i, T), possesses a pair of minima separated by a barrier 
for some range of temperatures [6]. Thus, we can classify models of the electroweak (EW) 
sector based on what physics is responsible for providing the requisite barrier in V^g (/i, T) . 
When calculated perturbatively, Kff (/i, T) is given by a sum of tree-level, quantum (loop), 
and thermal contributions. Thus, three model classes can be identified^ (see also Fig. 1): 

I. Thermally (BEC) Driven. A barrier arises due to thermal loop effects associated with 

bosonic zero modes. The effective potential acquires a term which ideally has the form 
— r(/i^)'^/^ where h is the Higgs condensate. Because the non-analyticity can be traced to 
the lowest energy mode of the Bose -Einstein distribution, this can be intuitively also called 
the Bose-Einstein condensation (BEC) driven scenario. The non-analytic term competes 
with the h? and h'^ terms in the scalar potential to generate a barrier. 

II. Tree-Level Driven. A barrier arises due to a competition between terms in the effective 
potential which are already present at tree-level. This model class can be further subdivided. 

IIA. Renormalizable Operators. The barrier arises from the competition between renor- 
malizable operators. Since an effective h? operator and gauge invariance are required for this 
class of models, these models rely upon a scalar field (or fields) in addition to a single Higgs 
doublet acquiring a nonzero expectation value during the EWPT. 

IIB. Non-Renormalizable Operators. If non-renormalizable operators involving the 
Higgs field (such as h^) are added to the scalai" potential, a barrier can arise as a result 
of their competition with the renormalizable terms. 

' However, it may be possible to weaken the tension between the Higgs mass measurement and the 
baryon asymmetry washout condition in non-standard cosmologies [-1, 5]. 

^ We do not claim that this classification scheme is exhaustive. For instance, models which rely on non- 
perturbative effects cannot be classified in this way (see, e.g., [7]). However, this classification does cover 
most perturbative models in the literature known to us. 
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I. Thermally (BEC) Driven IIA. Tree-Level (Ren.) Driven 




Higgs Field [ h ] Higgs Field [ h ] 



Figure 1. The four methods of obtaining a strongly first order phase transition by inducing a 
barrier in the thermal effective potential, which are discussed in this paper. The framed expressions 
indicate which term is responsible for the rise or fall of Ves- 



III. Loop Driven. Some h loop corrections may generate qualitatively important non-polynomial 
field dependence and aid in generating the barrier. For example, Ref. [ ] utilizes the quartic 
correction of the form In /i^, which can compete with the naively unstable — /i^ term to 
generate a barrier. 



In addition to a barrier in Vcs, successful EWBG requires the EW sphaleron process 
to be out of equilibrium in the broken phase to ensure that the baryon asymmetry is not 
washed out. This condition is expressed as a bound on the EWPT order parameter [?] 

^>1.3 (1.1) 

c 

T 

where {H)j, = (O , v{T)/ \/2) is the expectation value of the Higgs at temperature T, 
and Tc is the temperature at which the phase transition takes place (i.e., the symmetric and 
broken phases have degenerate free energy densities). We say that phase transitions which 
satisfy Eq. (1.1) are "strongly" first order phase transitions (SFOPT).^ 



^ The exact numerical value of the right hand side of Eq. (1.1) is mildly model dependent, but in all 
known cases, it is a number close to unity. 
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Thus, we will study the EWPT in the context of each model class by first parameter- 
izing the approximate thermal effective potential Vcs appropriate for each model class and 
then investigating what parametric limit will yield v(Tc)/Tc ^ 1. We can then ask what 
underlying physics would give rise to such an "optimal limit," what does the associated 
phenomenology look like, and what is the impact of collider constraints, assuming that 
the last Higgs-sector related phase transition is the electroweak symmetry breaking phase 
transition (i.e. there were no phase transitions that jumped from one electroweak symmetry 
breaking vacuum to another electroweak symmetry breaking vacuum). One of the conclu- 
sions of our study is that the optimal limits frequently correspond to enhanced symmetry 
points in the theory space. This makes the optimal limits straightforward to identify and 
associates them with a distinctive phenomenology which is constrained by recent LHC 
data. 

For example, one of the EWBG parametric regions most cleanly ruled out by the 125 
GeV Higgs is the enhanced continuous symmetry point parametric region (as opposed to 
the enhanced discrete symmetry point) which is a subset of Class IIA (Tree-Level Renor- 
malizable Operator Driven) models. As emphasized in [10], strong first order phase transi- 
tions can generically be found near parametric regions surrounding an enhanced symmetry 
point where the symmetry transformations mix Higgs and another field degree of freedom. 
One subset of enhanced symmetries is based on continuous symmetries (or the parametric 
limit in which the discrete symmetry enlarges into a continuous symmetry). One way to 
understand how the Higgs data rules out this subset is to note that the Nambu-Goldstone 
bosons associated with the spontaneously broken continuous symmetries have couplings 
to Higgs determined by the kinetic part of the action, and this coupling-induced decay rate 
is unsuppressed when the Higgs mass is of the order ofv = 246 GeV. Hence, the Higgs 
decay to the Nambu-Goldstone bosons exceed the experimental limits on exotic decays of 
the Higgs. 

The tension that we present in most of the categorization points to the enhanced dis- 
crete symmetry point [ ] being the parametric space marker having intuitively the largest 
set of model building possibilities for electroweak baryogenesis. 

In addition to constraints coming from the SM-likeness of the Higgs, it is also inter- 
esting to consider the "anomalies" which may point to beyond the Standard Model (BSM) 
physics. One of the most promising anomalies observed at the LHC is an excess of events 
in the loop-induced diphoton decay channel of the Higgs. If the excess can be attributed 
to the presence of a BSM scalar field running in the loop, then we utilize our classification 
to argue that there is a general tension with electroweak baryogenesis if this scalar field is 
also responsible for driving a SFOPT. 

The order of presentation is as follows. We begin with a review of the collider data 
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relevant for analysis. In Sec. 3, we present our classification and a general discussion of 
how the current data affects the models in the classification. We also give explicit model 
examples fitting into the proposed classification. In Sec. 4, we discuss the impact of the 
diphoton excess anomaly on each of the classes. We then close the paper with a conclusion. 

2 Collider Data and Interpretation 

Since models of the EW sector with strongly first order EW phase transitions tend to rely 
on a large coupling between the Higgs and light scalar fields, it is important to review the 
relevant constraints here. The Tevatron signal and ATLAS/CMS discovery confirm the 
existence of a bosonic particle with an approximate mass of 125 GeV [1, 3]. The available 
statistics suggest that the decays of this boson are consistent with the SM predictions in the 
channel 6H11, 12] as well as ZZ U and WW iviv [2, 3]. In the diphoton decay 
channel, both ATLAS and CMS observe an excess of events above the SM prediction 
[2, 3]. 

Spectrum: It is well known that in models such as the SM and the Minimal Super- 
symmetric Standard Model (MSSM), even the LEP Higgs mass bound imposes strong 
constraints on the viability of EWBG. Of course, these constraints have already ruled out 
EWBG in the SM [ ]. The measurement of a Higgs mass of 125 GeV further severely 
restricts the allowed MSSM parameter space [14, 15], although EWBG in MSSM is still 
viable with more judicious choices of parameters [16]. 

Exotic (e.g., Invisible) Decay: The discovery of a SM-like Higgs at the LHC is at 
tension with a large branching fraction in exotic channels. For instance, if the Higgs had 
a large branching fraction to invisibles, BRinv = BR(/i — > inv), this would suppress the 
branching fraction in all visible channels, and it would have been more difficult to find the 
Higgs at the LHG* [18]. A number of groups have investigated this possibility by assuming 
that the production cross section is the same as for a 125 GeV SM Higgs, but allowing for 
BRinv to vary in fitting the data. They obtain upper bounds on the branching fraction 
to invisibles in the range BRinv < 0.30 - 0.75 at 95% GL [i9-24]. Although this may 
not seem overly restrictive, we will see that in the phase transition model classes which 
allow invisible decay, this is naturally the dominant decay channel. Furthermore, the LHC 
expects to resolve the issue of invisible decay with increased data. It is estimated that 
with at 20fb~^ integrated luminosity, the LHC should detect or exclude invisible decay 

Assuming that the new physics does not enhance the Higgs production cross section, i.e., we assume 
(t(pp — > /i) = a^^iji-p — ^ K). However, even in the MSSM where new physics both allows invisible decay 
and enhances the Higgs production cross section, one finds that invisible decay is at tension with the data 
[17]. 
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for BRinv > 0.4 at 95% CL [25], and at 30 fb" , ATLAS should detect or exclude invisible 
decay for BRjnv > 0.24 at 5a ['' ]. Partially invisible final states resulting from exotic 
cascade decays are more difficult to constrain, but branching fraction bounds on the order 
of 10% may be obtained at the LHC with 1000 fb"^ [26]. 

Mixing with Hidden Sector: Just as with the case of invisible decay, the ATLAS/CMS 
data strongly constraints the scenario in which the Higgs is allowed to mix with a hidden 
sector scalar field or fields, which are singlets under the SM gauge group. For the sake of 
discussion, we will suppose that only one singlet scalar field is mixing with the SM Higgs. 
The impact of this mixing on the phenomenology depends on the relative mass scales, of 
the Higgs-like scalar atmn^ 125 GeV and the singlet-like scalar with mass mhid- Let 9 
be the angle between the Higgs-like mass eigenstate and the Higgs gauge eigenstate. The 
relevant constraints are: 

1. Light Higgs search at LEP. The existence of a light singlet-like resonance (i.e., 
rrihid ^ rriH = 125 GeV) is constrained by Higgs searches at LEP. In order for 
the singlet-like particle to have evaded detection, its coupling to the SM must be sup- 
pressed. This places an upper bound on 9, which which becomes more stringent as 
rrihid is deceased below the LEP Higgs search bound of 114.4 GeV. For instance, for 
mhid = 20 GeV one needs cos^ 9 > 0.99 at 95% CL [27]. 

2. Heavy Higgs search at LHC. Similarly, if the singlet-like resonance is heavier (i.e., 
rrihid ^ m-H = 125 GeV), there is an upper bound on 9 coming from the requirement 
that the heavy singlet-like scalar evades detection at the LHC. Again, this is a function 
of the singlet-like scalar's mass. For instance, if mhid = 200 GeV one needs cos^ 9 > 
0.60 at 95% CL to avoid detection [28, 29]. 

3. LHC Higgs Detection. Assuming that the Higgs-like resonance is lighter (i.e., 
rriH = 125 GeV <^ mhid)> then the consequence of mixing is a universal suppres- 
sion of all Higgs production processes by a factor of cos^ 9. Large mixing would have 
made discovery more difficult. Thus, the LHC's signal at 125 GeV places an upper 
bound on 9 which may be expressed as cos^ 9 > 0.77 at 90% CL [24, 30]. (See also 
[19, 28, 31-33]). 

Taken together, these constraints imply that the large mixing scenario (e.g., cos^ 9 = 0.5) 
is strongly disfavored. 



3 Electroweak Phase Transition IModel Classes 

In this section, we will enumerate the phase transition model classes, identify the paramet- 
ric limits which are optimal for SFOPT by maximizing the washout criterion Eq. (1.1), 
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and discuss phenomenological constraints that arise in those limits. To connect to phe- 
nomenological constraints, we will make a simplifying assumption that the electroweak 
symmetry breaking is the last Higgs-sector related phase transition (i.e. no transition from 
one electroweak non-symmetric vacuum to another electroweak non-symmetric vacuum).^ 
As we discuss further below, the optimal limits for SFOPT often correspond to enhanced 
symmetry points of the theory at which the symmetry group is extended to include an ad- 
ditional continuous or discrete symmetry. For the sake of brevity, we will not dwell on the 
details of the phase transition calculation. We refer the interested reader to the review [6]. 

3.1 Class I: Thermally (BEC) Driven 

In models such as the SM and the MSSM, the barrier in the thermal effective potential 
arises from thermal loo^ effects, which emerge in the following way. The Higgs conden- 
sate {H) = (O ,h/ \/2j modifies the dispersion relation of particles in the plasma causing 
them to acquire an effective temperature and field-dependent mass Tnlg{h, T) = m'^(h) + 
n(T). Here, 11 is a temperature-dependent self-energy correction (known as "daisy resum- 
mation," see e.g., [ ]) and fri^ih) can be obtained by replacing the zero temperature VEV v 
with h in the standard expression for the field's mass (see e.g., [ ]). Bosonic fields induce a 
contribution to the thermal effective potential of the form \4fr 3 {—T/\2-n) (m^g {h, T)) 
in the high-temperature limit. The non-analyticity of this term at m^g = can be traced to 
the non-analyticity of the Bose-Einstein distribution function at zero energy. Hence, this 
thermal "BEC term" driven SFOPT defines our "Class I" model class. 

To achieve a barrier in V^g near the phase transition temperature Tc, we want to have 
(m^g(/i, Tc))'^/^ ~ such that there may be a competition between this term and the h? 
and terms of the Higgs potential. Supposing that fh^ih) can be written as fh?{h) = 
ah"^ + (3, the effective mass will have the desired scaling if we tune (3 — n(Tc) ^ af (Tc)^. 
A general phenomenological consequence of this tuning is that the scalar bosons today 
will be light, since their mass squared is m'^{v) ^ av'^ — n(Tc). Note that increasing the 
interaction of the h field to make a large naturally drives up 11 (Tc) which in turn drives 
m^{v) lighter. The need for this tuning using (3 is well-established in the MSSM [35], 
where light right-handed stops provide this (m^g(/i, Tc))^^"^ term. Phenomenologically, the 
light stops tend to enhance Higgs production by gluon fusion and reduce Higgs diphoton 
decay. Because of this, the LHC has already placed strong constraints on EWBG in the 
MSSM [14, 15], and it has begun to push the model into a comer that will be probed by 
the high luminosity LHC [ 1 6] . 

^With a sufficiently large number of broken vacua jumps in between, the electroweak symmetry breaking 
vacuum properties measured at colliders today can be decoupled from the vacuum properties associated with 
the first electroweak symmetry breaking phase transition. 
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Near the temperature of the phase transition, the effective potential may be approxi- 
mated as 



VesiK T) ^ \ (V + cT^) - ^ih^' + (3.1) 

in the high temperature expansion. Note that a factor of l/127r has been included in the 
parameterization to reflect the natural thermal loop suppression of this coefficient. A po- 
tential of this form is illustrated in Fig. 1. The parameters =m|^/2andA = m|^/(2f^) 
are related to the Higgs mass tuh and VEV v.^. The dimensionless parameters c and e 
quantify the coupling between the Higgs condensate and the relativistic particles in the 
plasma. In particular, c depends on couplings between h and light (m < T) bosons and 
fermions, whereas e only depends upon couplings between h and light bosons. Schemati- 
cally, 

e ~ (degrees of freedom) x (coupling to Higgs) ^''^ . (3.2) 

light bosonic fields 

The contribution from heavy fields (m > T) are Boltzmann suppressed, and the 0(T^ exp [—m/T]) 

terms are dropped. Some examples of models that fall into this class are shown in Table 

1^. 



A standard calculation (see, e.g., [ ]), yields the EW order parameter 

v(Tc) e 



Tc QttX 

There are two "optimal" limits in which we can obtain v{Tc) /Tc ^ 1. 



(3.3) 



e ^ A To reach the limit of large e, the Higgs must have a large coupling with many 
light bosonic degrees of freedom. Indeed, the presence of Gtt in the denominator of 
Eq. (3.3) (which comes from the thermal loop expansion) makes satisfying Eq. (1.1) 
very challenging if A ~ 0(1). There are various phenomenological constraints on 
this limit. First, since e is a sum of dimensionless coupling constants (see, e.g.. 
Table 1), it is bounded from above by the perturbative unitarity constraint. Second, 
heavy bosonic fields will become Boltzmann suppressed and cannot contribute to 
e. However, the same interactions which allow light bosonic fields to contribute 
to e also provide a mass to those fields after EWSB. Thus, increasing the coupling 
constants that enter e, will eventually cause the bosons to become heavy and their 

''More generally, h need not be the Higgs and mn need not be the 125 GeV Higgs mass. 

^In models such as the MSSM and "Colored Scalar" model, the light scalars that provide the BEC 
term are colored (e.g., stops in the MSSM). Higher order contributions, which scale like —g'^(N^ — 
l)T^m? log mi/T, tend to lower the PT temperature and strengthen the PT 
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Model 
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SM r , 1 




6m? +6m?,,+3m?, + 4m|, 
rem = w z 1 H 


6m,?,, +3r?i?. 
Pour = „ ^ 


MSSM 

[-33] 




, 6m2 A2 \ 
CSM + 12^1 (^1 - 4 j 


, 6m3 Al 


Colored 
Scalar [15] 


Ml 
Q\H 


X 

-4 
2 


f + 
+ 


„ 1 6 Q 
CSM + 24 2" 


esM + 6(f)^/^ 


Singlet 
Scalar [37] 


M^\S\ 

2e \H\ 


' + 

+ 


+ f,e 


esM + gsC 


Singlet 

Majoron 

[38] 


A. 1^1' + 
Xhs\Hf\Sf+ 
lyiSPiUi+h.c. 







Table 1. Examples of models in the Thermally (BEC) Driven class. The expressions for e are 
calculated in the limit that the field independent contributions to mlg{h, T) are negligible (e.g., the 
thermal mass tuning has been performed). Here, the symbol At is At = At — n/ tan 13 and gg is 
the number of real scalar singlet degrees of freedom coupling to the Higgs. 

contributions to e will become Boltzmann suppressed.^ (One can however increase 
e up to the perturbativity bound by increasing the number of degrees of freedom that 
contribute to the thermal loop instead of increasing the coupling constant. However, 
in that case, one may need to arrange fermion loops to cancel radiative corrections to 
A.) Finally, as e is increased, interactions between the Higgs and other bosonic fields 
are made stronger. Thus, there may be loop-suppressed, but nevertheless significant, 
modifications to Higgs production and / or decay. For example, if the bosons carry 
color, then they can significantly enhance Higgs production by gluon fusion [15, 39]. 
We will revisit this constraint in the context of Higgs diphoton decay in Sec. 4. 

A — >• In the context of the SM, this limit is obviously forbidden in light of the re- 
lationship A = mjj/2v'^ and the fact that uih is now a measured quantity. However, 
in an effort to keep our model classification scheme as general as possible, we will 
consider the scenario in which the field h that appears in Eq. (3.1) is not the SM 
Higgs condensate. Instead, it may represent a parameterization of some non-trivial 
trajectory through an extended scalar field space connecting the EW-preserving vac- 
uum h = with the EW-broken vacuum h = v. Then, the limit A — )■ implies 

^Note that the examples we have chosen in Table 1 do not include Z2 breaking cubic couplings since 
those would naturally have strong phase Uansition possibilities driven by tree level terms. 
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the spectrum contains a light scalar. If the scalar carries SM quantum numbers, then 
direct search constraints are severe unless the scalar can be hidden in a large SM 
background. If the Higgs decay channel is open, this limit may be at tension with 
constrains on Higgs exotic decay. If the scalar is a SM singlet, then constraints on 
hidden sector mixing may also apply. Moreover, vacuum stability considerations 
limit the range of the EFT (see, e.g., [ ] and references therein). 



To illustrate how these limits and constraints arise in a concrete model, we extend the 
SM by a color triplet scalar field X (see [15]): 



£ = £sM + {d^Xy 



MlX*X 



K 
6" 



{X*Xf + QH^HX*X 



(3.4) 



Note that the quantum numbers of X only allow it to couple to the EW sector via the so- 
called "Higgs portal" operator H^X^ with coefficient Q. The effective mass of the scalar 
X is given by m\{h,T) = Mj- + {Q/2)h^ + IlxiT) where Ilx{T) = {K + Q)Ty24. 
Thus, the BEC term is given by 



AV;ff(/i,T) = -6 



T 

12^ 



Mi 



nx(T)j 



\ 3/2 



(3.5) 



where the factor 6 is the number of internal degrees of freedom for the complex, colored 
X field. As discussed above, we must tune ~ —Ilx{Tc)- Thus, Eq. (3.5) takes the 
form of Eq. (3.1) with e = csm + 6(Q/2)3/2. 



We would like to understand what constraints arise as we go to the SFOPT limit e ^ 
A. This limit is reached by taking Q ^ X^^^. First, we verify that the phase transition is 
strongly first order by calculating v{Tc)/Tc as a function of Q. We fix mu = 125 GeV, 
K = 0.1, and require M|- = -Ilx{Tc) where 

T = [l - — (—] 

The numerical calculation is performed in the standard way (see, e.g., [6]) using the full 
one-loop, daisy improved thermal effective potential. As shown in Fig. 2, the EW phase 
transition becomes strongly first order for sufficiently large values ofQ > 1.7. Second, we 
note that perturbativity up to 100 TeV requires Q < 2 at the weak scale [ ]. Third, as we 
discussed above, Boltzmann suppression of heavy X bosons prevents us from obtaining 
SFOPT for arbitrarily large Q. We can estimate an upper bound on Q by requiring the 
X bosons to be light at the temperature of the phase transition: Tc > mxiv{Tc),Tc) ~ 
^yQ/2v{Tc) translates into Q < 2(Tc/f (T^))^ ^ 2. The numerical calculation confirms 
this estimate and explains the discrepancy between the numerical an analytic calculations 
ofv{T,)/T,. 



(3.6) 
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Figure 2. A plot of the EW order parameter v(Tc)/Tc calculated analytically (black, dashed) and 
numerically (black, solid), as discussed in the text, as well as modifications to Higgs production by 
gluon fusion (red) and Higgs decay to two photons (purple). The numeric calculation of v{Tc)/Tc 
falls short of the analytic estimate due to the Boltzmann suppression effect discussed in the text. 
The analytic expression suggests that SFOPT are obtained for Q > 1.2, but numerical calculation 
reveals that SFOPT are only found for Q ^ 1-7. Hence, there is a narrow window 1.7 < Q < 2.0 
where the perturbative calculation is valid and the EWPT is strongly first order. In this region, 
when X is a electrically charged color triplet, the phenomenology consists of an enhanced rate of 
gg ^ H and a reduced i7 — )• 77 rate. 



As noted in [' the addition of a color triplet X field to the SM is motivated by a 
desire to use the 6 real scalar degrees of freedom to enhance the strength of the EWPT. 
Such a scenario has a natural connection with collider physics, via the X field's contribu- 
tion to the rate of Higgs production by gluon fusion (see also [41]). Furthermore, if the 
X field is electrically charged like the stops, then it can also affect the diphoton decay 
rate of the Higgs. To illuminate this point, we also show on Fig. 2 the modifications to 
gluon fusion and diphoton decay, which are calculated following [42] and using gx = 2/3 
for the electric charge of the X field. As shown in Fig. 2, in the SFOPT window, gluon 
fusion is enhanced and i7 — )► 77 decay is suppressed by an 0(1) factor with respect to 
the SM rates. This phenomenological signature is marginally disfavored by the recent 
ATLAS/CMS data, but additional statistics will be required to justify a strong statement. 

In summary, even if the X fields are SM singlets and thereby able to evade collider 
restrictions. Fig. 2 illustrates a generic strong tension from phase transition and theoretical 
considerations alone. The tension is ultimately tied to the difficulty of overcoming the nat- 
ural 1 /Gtt suppression appearing in Eq. (3.3) while remaining in the perturbative regime. 
The most promising way of overcoming this natural suppression is to have a small ef- 
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Model 




xSM[ J, 44] 


I {dSf - [^S^ + + 'fS^ + fWHS^ + fH^HS^] 


ZsxSM [10] 


I (dSY - [1^2 + + fWHS^] 


Model 


AW 


NMSSM [4j] 


XH1H2N - f iV^ 


nMSSM ["-^ 


XHiH2S+'^S 


/iz/MSSM [ ] 





Table 2. Examples of models that fall into Class IIA. For the non-SUSY models, corrections to the 
SM Lagrangian are shown, whereas for the SUSY models only the superpotential corrections are 
given. 

fective quartic coupling which is model-dependently constrained by collider observations 
since it typically signals light particle states which have not been observed. 

3.2 Class IIA: Tree-Level (Renormalizable Operators) Driven 

We saw in the previous section that the Thermally (BEC) Driven models are strongly 
constrained, ultimately because of their reliance on the BEC term and its thermal loop 
suppression factor of l/Gvr. Our next class of models which we call "Class 11" relies 
instead on tree-level interactions of the Higgs to provide the barrier for the SFOPT. For 
renormalizable models, these tree-level operators are cubic in the fields (with respect to 
a particular field origin associated with the EWPT). Then, gauge invariance require that 
there be at least one scalar in addition to the SM Higgs that acquires an expectation value 
during the EWPT. For non-renormalizable models, a barrier may be obtained without any 
odd powered monomial terms, and therefore we will further subdivide the tree-level model 
class into two subclasses ("Class IIA" and "IIB"). We will find perhaps the most clean 
nontrivial result of this paper that a particular comer of the Class II model class is ruled 
out due to the current Higgs data. 

First, we consider the class of models (which we call "Class IIA") in which the barrier 
in VcQ arises from renormalizable tree-level interactions between the Higgs and new scalar 
fields. Thus, the term in V^s that provides the barrier is necessarily cubic, and derived 
from dimension 3 or 4 scalar interactions in the Lagrangian. This naturally does not suffer 
from the l/Gvr thermal loop factor handicap as in Class I models. As we remarked above, 
at least a single scalar degree of freedom in addition to the SM Higgs must participate in 
the phase transition. We thus parameterize the additional scalar field(s) as S. The number 
of degrees of freedom associated with S, its quantum numbers, and its interactions will 
be model-dependent. The information that is pertinent to our generic phase transition 
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analysis is that there exists a one-dimensional trajectory through the configuration space 
which interpolates between the EW- symmetric and EW-broken phases. The effective 
potential along this trajectory may be approximated as 

where we have only included the leading high temperature dependence, since by definition 
of this model class, we are assuming that the temperature independent tree-level Lp^ is more 
important than the naturally suppressed TiyLp^Y^"^ /12'k term. The fact that this model class 
does not have to generate a cubic term dynamically and overcome the natural l/127r sup- 
pression gives this class a considerably larger model freedom than Class I. Some examples 
of models that fall into this class are shown in Table 2. 

The phase transition temperature is calculated from Eq. (3.7) to be 
and the EW order parameter is found to be 



V A / 1 \m? 



COS a. (3.9) 



Here, we have introduced a projection factor of cos a, since in general if will not corre- 
spond to the Higgs field. The optimal limits for enhancing v(Tc)/Tc are given by: 



c ^ A Since c represents a sum of coupling constants controlling interactions be- 
tween the Higgs and light particle in the plasma, one might try to take the limit 
c ^ A by increasing the size of these couplings or by increasing the number of 
degrees of freedom in the plasma. Although this limit is similar to the e ^ A case 
discussed for Class I, they differ significantly in that e only receives contributions 
from bosonic degrees of freedom (recall the name BEC driven), whereas c receives 
contributions from fermions as well. The Higgs self-coupling A is also renormalized 
by these same couplings that enhance c. Generally, it is not obvious that the limit 
used to increase c will not also increase A and thereby prevent one from reaching the 
c ^ A limit. For example, we can consider the contributions to c and A that arise 
from the Yukawa interaction with the top quark. The contributions scale with the 
Yukawa coupling ht and number of colors Nc like c ~ N^hf and A ~ —NJil yield- 
ing c/A ~ —\/hl. In this example, increasing the value of the Yukawa coupling 
will tend to decrease the ratio of c/A. One way to get around this result is to note 

^We can parameterize the one-dimensional trajectory with a field cp, as h = h{(p, T) and S — S{ip, T). 
In principle, the functions h and S can be determined by solving for the multi-field bounce solution. 
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that contributions to c are non-negative whereas contributions to A are positive for 
bosonic fields and negative for fermionic fields. If the underlying model possesses 
a symmetry relating bosonic and fermionic fields (such as SUSY) then it may be 
possible to take c large while keeping A small. If the light fields do not carry any 
SM quantum numbers, and if they are sufficiently light (m < m^/S) then \ is 
at tension with constraints on Higgs invisible decay. 

Xm^/2£^ — >■ 1 This is the limit in which Tc vanishes and the EW-symmetric and 
EW-broken vacua are degenerate. As noted in [ ], this degeneracy may arise as the 
result of a discrete symmetry relating the Higgs field with the other field(s) partic- 
ipating in the phase transition. We will refer to this limit as an enhanced (discrete) 
symmetry point (EdSP), which is illustrated in Fig. 3. As one approaches the EdSP, 
the EW-symmetric vacuum becomes metastable and increasingly degenerate with 
EW-broken vacuum. Without sufficient degeneracy breaking, tunneling out of the 
EW-symmetric vacuum may become suppressed to the point that tunneling occurs 
on a time scale that exceeds the age of the universe. That is, as one approaches 
the EdSP, it may be the case that the EWPT never occurs, even if the EW-broken 
vacuum is energetically favored. 

A — >■ We would like to take this limit while fixing Xw? /2£'^ such that Eq. (3.9) 
just scales like 1 / \/\. Moreover, if we also want to fix the VEV of the (p field 



then we see that we must let £^ oc A and oc A as A — 0. In this limit, the mass of 
the ip field 



also scales like A and goes to zero. Thus, there will be a light scalar field associated 
with the Lp field direction. The light scalar runs into two phenomenological con- 
straints. If ip represents a mixture of the Higgs with a hidden sector scalar field, then 
a light Higgs is excluded by searches at LEP and at tension with the LHC Higgs dis- 
covery. On the other hand, even if there is no mixing, provided that the light scalar 
is mostly a SM singlet, then this limit runs into constraints on Higgs invisible decay 
imposed by the LHC Higgs discovery. We will discuss this scenario further in an 
example below. 

It is important to note that as we take this limit in which m?, 8, and A approach zero, 
the effective potential develops a shift symmetry. Thus, we can identify the A — i- 
limit with a enhanced symmetry point of the theory at which a continuous symmetry 




(3.10) 




(3.11) 
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Figure 3. A illustration of the behavior of Vcs as the limits Xm?/2£'^ — )• 1 (left) and A — >• (right) 
are taken. The former leads to an EdSP whereas the latter leads to an EcSP. 



emerges. We will refer to this parametric limit as a enhanced (continuous) symmetry 
point (EcSP), which is illustrated in Fig. 3. 

In order to demonstrate how these limits and constraints may be realized in a concrete 
model, we consider the Z2XSM [ ]. This model extends the SM by a real scalar field S 
which is a singlet under the SM gauge group, but which respects a Z2 discrete symmetry 
that takes S — )■ —S. The most general, renormalizable Lagrangian consistent with the SM 
gauge group and Z2 is given by'" 



■^Z2xSM — -^SM + 2 i^fj."^) (d^'S) — 



+ + "^H^HS' 



(3.12) 



where Csm is the Lagrangian of the SM. We assume that S does not acquire a VEV. Thus 
the Z2 is unbroken, thereby ensuring stability of S and preventing mixing with the Higgs. 
Although S does not have a VEV, we will allow it to obtain a non-zero expectation value at 
finite temperature so that it may participate in the EWPT and render it strongly first order. 



With this Lagrangian, we can calculate the effective potential as a function of both the 
Higgs condensate (H) = (O , h/\/2j and singlet condensate {S) = s. Working to the 
same level of approximation as in Eq. (3.7), we neglect the loop-suppressed contributions 
and include only the leading thermal contributions to obtain 



Vesih,s,T) 



: 4 



-62 + C.T2 2 



4„4 , «2^2^2^ 



+ — + 

4 4 



(3.13) 



The thermal mass terms Ch and Cg ensure symmetry restoration at sufficiently high 
temperature. 



'"since the one-loop phase transition analysis does not depend upon the quantum numbers of S, the 
analysis here will also apply to the more general case of a non-singlet S coupled via the "Higgs portal" 



operator W HS* S. Such a scenario is discussed in Sec. 4. 
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In light of the general analysis of the preceding subsections, we are motivated to seek 
out enhanced symmetry points. In the following discussion, we will identify the EcSP, 
justify the claim that SFOPT are found in its vicinity, determine the phenomenology in 
this limit, and assess the impact of collider constraints. We will then repeat the analysis 
for a neighborhood of the EdSP. 

The parameters of the Z2XSM are the SM gauge (gi) and Yukawa couplings (yi), the 
Higgs sector parameters (/i^ and A), the singlet sector parameters (62 and 64), and the 
"Higgs portal" coupling (02). The symmetry group of the Z2XSM Lagrangian is Gsm x 
Z2 where Gsm is the gauge group of the SM. For a particular choice of parameters, the 
symmetry group enlarges to incorporate an additional continuous symmetry. We find this 
EcSP by requiring 

EcSP : {&2 = /i^ , &4 = A , 02 = 2A} and {gi = , Vi = 0} , (3.14) 

where A = ni\j/{2v'^) and jj? = are not constrained by the symmetry, but are 

restricted by measurements of the Higgs mass and VEV. At the EcSP, the Lagrangian can 
be written as 



SM 



EcSP Z 



(3.15) 



up to kinetic terms for the other SM fields. By virtue of the EW symmetry, this Lagrangian 
is invariant under an SO (4) symmetry which acts on the components of H. However, by 
virtue of the EcSP, this symmetry is enlarged to an S0(5) group which rotates among 
the components of H and S. The symmetry ensures that Cs = Ch = cq and the effective 
potential may be written as 

s^T) =\ (-/.^ + CO T^) {h' + s') + ] {h' + s')\ (3.16) 
EcSP 2 ^ ^ ^ ^ 4 ^ ^ 

Evidently the restriction to vanishing gauge and Yukawa couplings is unphysical, and once 
these couplings are turned on, radiative corrections to Vcs will break the SO (s) symmetry 
back down to SO (4). However, the symmetry breaking terms will carry a loop suppres- 
sion factor of l/lGvr^ and can be neglected at this level of approximation. On the other 
hand, contributions to the thermal masses are not loop-suppressed and will generically 
induce Ch 7^ Cg. Therefore, in the following discussion we will neglect loop suppressed 
corrections to Vcs, but we will treat Ch and Cg as independent parameters. 

We will see that there are SFOPT in a neighborhood of Eq. (3.14), but first it is inter- 
esting to remark that the pattern of symmetry breaking is controlled by the symmetry that 
arises at this EcSP. Provided that /i^ > 0, the continuous symmetry will be spontaneously 
broken. The resulting Nambu-Goldstone boson is associated with a flat direction in the 
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potential connecting \H\ = vl^2 with S = v. Thus, we anticipate that we will find phase 
transitions that occur in two steps: first S acquires an expectation value breaking the Z2, 
and second the expectation value of S returns to zero as H acquires an expectation value 
breaking the EW symmetry. 

We can proceed to perturb away from the EcSP by writing the parameters as 

62 = /i^ (1 + CbJ , 64 = A(l + e64), and 02 = 2A (1 + CaJ • (3.17) 

What sort of perturbations will yield SFOPT? At the EcSP, the EW-broken and EW- 
symmetric vacua are degenerate, and if Ch = Cg then the thermal corrections will maintain 
that degeneracy. As we perturb away from the EcSP looking for SFOPT, we will need to 
ensure that degeneracy breaking causes the EW-broken vacuum to be energetically favored 
and also ensure that the breaking of Ch 7^ Cs causes the EW-symmetric vacuum (in which 
Z2 is broken) to become (free-)energetically favored above some temperature. Keeping 
this picture in mind, we can proceed to calculate the phase transition parameters. In this 
neighborhood of the EcSP, the phase transition temperature and EW order parameter are 
given by 

= ^J!l^ ^e,,-2et,{l + 0(e,„ e,j) (3.18) 

^ = 2V^r^— ^ ^ ^ (1 + 0{eb„ e^j) . (3.19) 

See also Fig. 4. As we anticipated, Tc is arbitrarily small and v{Tc)/Tc is arbitrarily large 
for arbitrarily small perturbations away from the EcSP (e;,^ — 2eb2 ^ The particular 
combination of parameters e^^ — 2eb2 appears, because it controls the degree of degener- 
acy breaking between the EW-symmetric and EW-broken vacua. We can verify this by 
calculating 

Kff {0,Vs,T)-Ves {v,0,T) 

^ (e,, - 26,J -^{Ck- Cs) T'] (1 + 0(e,2 ~ e,, ~ T^)) , (3.20) 



where Vg = ^/b2/b4 is the expectation value of s in the EW-symmetric vacuum and v = 
a/ yU^/A = 246 GeV is the Higgs VEV. Thus if e^^ — 26^^ = 0, the two vacua are degenerate 
at T = 0. If ef,4 — 2efc2 > 0, the broken vacuum is energetically favored and the PT occurs 
at the temperature Tc given by Eq. (3.18), but if e^^ — 26^2 < 0, the symmetric vacuum 
is energetically favored and the PT does not occur. From this discussion, and particularly 
Eq. (3.19), we conclude that SFOPT are found in the neighborhood of the EcSP, but 
additionally the EcSP demarcates a boundary between physical models (e^^ — 2efc2 > 0) 
in which EWSB occurs and unphysical models (e^^ — 26^2 < 0) in which EWSB does not 
take place. The singular factor of 1/ y/ch — Cg in Eq. (3.18) can also be understood in light 
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Figure 4. SFOPT correlated with large invisible decay in a neighborhood of the EcSP. The dashed 
lines corresponds to values of the singlet mass ms- Left. The EW order parameter v{Tc)/Tc, for 
which Eq. (3.19) is the leading order expression. Right. The branching fraction of Higgs to an 
invisible S-pair BRjnv, for which the width Eq. (3.22) is the leading order expression (see also 
[18]). 



of Eq. (3.20). If Ch = Cg, then thermal corrections lift the EW-broken and EW-symmetric 
phases together maintaining their degeneracy. One needs Ch > Cg to ensure that Vcs at the 
EW-broken phase (free energy density) is lifted more greatly with increasing temperature 
than the EW-symmetric phase. Conversely, if Ch < Cg then the EW-symmetric phase in 
which the Z2 is broken never becomes (free-)energetically favored. 

We can begin to investigate the phenomenology near the EcSP by calculating the mass 
of the singlet scalar field. The tree-level relationship can be read off of the Lagrangian 
Eq. (3.15), which gives 

2 t , «2 2 EcSP mjj /o 1 N 

ms = -&2 + ^ ~Y~ ^^"^ ~ ^''2) • (3-21) 

Since this scalar field corresponds to the pseudo-Nambu-Goldstone boson of the sponta- 
neously broken continuous symmetry, we are not surprised to find that it is light when 
deviations away from the EcSP are small e^^ — ^ 1. Note that at the EcSP, the Nambu- 
Goldstone boson kinetic term will couple to the Higgs through a dimension 5 operator 
with a coupling strength of order f /v where / is a group theory factor typically of order 
unity. This generically leads to a large Higgs invisible'' decay width. In the toy model at 



" More generally, Z2-violating couplings between the hidden sector and the SM may allow S to decay 
back into SM particles. In that case, the same constraints apply to the unobserved exotic decays. 
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hand, the decay width is 

T{H ^ SS) ^ + (6,, + + 0{el,el)) . (3.22) 

See also Fig. 4. Since S only couples to the SM via the Higgs, the width for Higgs 
decay into SM fields, T(H ^ SM), is only affected by its coupling to S at the multi-loop 
level. Thus, we can approximate T(H SM) by the SM Higgs total width, which is 
~ 5 MeV for ttih ~ 125 GeV [ ]. We find that the invisible branching ratio is 

where we also neglect kinematically suppressed 3-body (and greater) final states. Such a 
large invisible decay greatly exceeds the 95% CL limits set by analyses of the LHC and 
Tevatron Higgs data, which were discussed in Sec. 2. Thus, the tension which we had 
discussed between the EcSP limit and invisible decay is illustrated in a concrete setting. 

We can attempt to evade the collider constraints on Higgs invisible decay by suppress- 
ing the channel H — )■ SS. This can be accomplished by moving away from the EcSP. 
In the following, we will discuss two ways of deviating away from the EcSP while main- 
taining a SFOPT. The first way will be to reach an enhanced discrete symmetry point in 
the parameter space such that the mass of the would-be Nambu-Goldstone boson is lifted 
above the threshold for two-body decay of the Higgs. The existence of a remnant symme- 
try is what guarantees the SFOPT in this first deformation [ ] . The second way will be to 
approach a free theory limit for the would-be Nambu-Goldstone boson while maintaining 
a symmetry of the potential at the tree level. In this second deformation, the kinetic term of 
the would-be Nambu-Goldstone breaks the symmetry of the potential, but such breaking 
is mild enough to ensure a SFOPT [ ] . 

Let us consider the first deformation. Specifically, we perturb away from the EcSP 
such that the continuous symmetry is broken to its discrete subgroup §2 which exchanges 
y/2H o S. The EdSP is given by 

EdSP : {62 = yu' , h = \} and {g^ = , = 0} . (3.24) 

Since /i^. A, and 02 are free to vary, the EdSP represents a 3 -dimensional submanifold of 
the full Z2XSM parameter space in contrast to the 2-dimensional submanifold correspond- 
ing to EcSP. As before we can consider perturbations away from the EdSP parametrized 
as 

62 = /i^ (1 + CfcJ and b4 = X{l + ebJ. (3.25) 

Because the singlet is no longer the Nambu-Goldstone boson of the spontaneously broken 
continuous symmetry, its mass need not be small: 

2 . ^2 2 EdSP oxx A 2A \ 

ms = -b2 + Y' - 2A) (1 - ^^^e,, 1 . (3.26) 
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From this expression we can see how the variation of 02 affects the vacuum structure. 
For a2 = 2A we return to the EcSP and the singlet is the massless Nambu-Goldstone boson. 
For 02 < 2A, the singlet becomes tachyonic, signaling that the true vacuum of the theory is 
one in which the Z2 is spontaneously broken. This is an undesirable limit, because without 
the Z2 preventing the Higgs and singlet from mixing, we run into the collider Higgs mixing 
constraints, which were discussed in Sec. 2. For 02 > 2A, the vacuum preserves the Z2 and 
the singlet is massive. Provided that 02 > 3 A, the singlet mass ms > ^1^/2 will exceed the 
kinematic threshold and block the invisible decay H — )■ SS. Using mn ~ 125 GeV and 
A = mjj / {2v'^) , this bound is approximately 02 > 0.39. Moreover, since the expressions 
for the phase transition temperature and EW order parameter, Eq. (3.18) and Eq. (3.19), 
were independent of 02, we still expect to find SFOPT in this corner of parameter space 
near the EdSP. Thus a departure from the EcSP along the EdSP allows for SFOPT while 
avoiding Higgs invisible decay by kinematically blocking the H SS channel. 

A second deformation away from the EcSP while preserving SFOPT but avoiding 
Higgs invisible decay is obtained by moving towards a free theory (02 = 0) for the would- 
be Nambu-Goldstone field without making it heavy. To maintain the SFOPT, we must do 
this while preserving the degeneracy of the energy of the two vacua involved in the phase 
transition. In the previous discussion we saw that if we moved away from the EcSP along 
the direction of the EdSP, then taking 02 < 2A would lead to an undesirable change in 
the vacuum structure such that the Z2 becomes spontaneously broken and the Higgs and 
singlet are allowed to mix. Thus, we must find a different path that continuously connects 
the EcSP with 02 = but maintains the vacuum structure including the degeneracy. 

The correct path is given by the following parameter choice. 



EcSP: = , 64=(^j A|. (3.27) 

At the parameter point Eq. (3.27), the scalar sector Lagrangian can be written as 



Z2xSM 



EcbP Z, 

-/i^ [h^H + ^572) + A (^H^H + IJ5V2VI . (3.28) 



From this expression we see that the scalar potential is invariant under a continuous sym- 
metry transformation which rotates and dilates the fields H and S, but that the scalar 
kinetic terms are not invariant (unless 02 = 2A). Thus, Eq. (3.27) is not a true enhanced 
symmetry point of the Z2XSM. Radiative corrections will spoil the symmetry, and there- 
fore we do not expect the effective potential to respect this symmetry (even if we were 
to also set Qi = Ui = 0). Nevertheless, since in this class of models, the phase transition 
parameters are dominantly controlled by the structure of the tree-level scalar potential, we 
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Figure 5. The EW order parameter f (Tc) /Tc (blue) and invisible branching fraction BRinv (red), 
calculated as in Fig. 4 but at the EcSP' parameter point Eq. (3.29). As a2/2A is decreased below 
about 0.15, Higgs invisible decay becomes sufficiently suppressed to evade collider constraints 
which impose BRinv ^ 0.64. This threshold corresponds to 02 ~ 0.043. At the same time, the 
electroweak phase transition remains strongly first order v{Tc)/Tc > 1. 

expect that SFOPT may still be found in the vicinity of Eq. (3.27). However, it turns out 
that in breaking this continuous symmetry the radiative corrections split the degeneracy 
of the EW-symmetric and EW-broken vacua in such a way that the EW-broken vacuum 
becomes metastable, and consequently EWSB does not occur. To avoid this outcome, we 
must allow for a finite perturbation away from the EcSP parameter point. We consider 
instead the EcSP' defined to be 



where we will allow 02 to vary and keep = e^^ = — 1/2. 

Along the trajectory Eq. (3.29) we can take a2 while keeping A and /i finite. The 
singlet remains light m| = —62 + 02^^/2 = a2v'^eb2 and the invisible width is approxi- 
mately given by 



To bring the invisible branching fraction below BRinv < 0.64 (one of the weakest 95% 
CL limits [ ]) we need 02 < 0.043. Furthermore, since the expression for the EW order 
parameter Eq. (3.19) is independent of 02, we still expect to finds SFOPT in this limit. 
This can be verified by calculating the EW order parameter numerically, and the result is 
shown in Fig. 5. 




(3.29) 




(3.30) 
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In summary, within the Class IIA scenario, EcSP region of the SFOPT parametric 
region is cleanly ruled out by the current data disfavoring large Higgs branching to BSM 
states. What is clear from the two deformations away from EcSP in the context of a simple 
BSM model is that EdSP does not require a small dimensionless parameter while EcSP' 
requires a tiny dimensionless coupling which begs for an explanation. In that sense, EdSP 
more naturally accommodates both a SFOPT and the current Higgs data. 

3.3 Class IIB: Tree-Level (Non-Renormalizable Operators) Driven 

The second way of obtaining a SFOPT using only tree-level operators is to employ non- 
renormalizable terms in the potential. If the scale of new physics A is not much larger than 
the EW scale, then the leading correction to the scalar potential, (WH)^, may dramatically 
change the nature of the EWPT.'^ In this scenario, the effective potential may be written 
as 



Since typically v{Tc) < v, the 0{h^/ A'^) terms can be neglected provided that A > t>. By 
minimizing the potential, the parameters /i^ and A may be exchanged for the Higgs VEV v 
and mass ttih- These relationships are given by 



where we have introduced Amax = y/Sv'^/fnn, the meaning of which will become clear 
shortly. Since we are interested in the limit that will yield a barrier in the effective potential, 
we will focus on the case of a low-scale cutoff such that /i^ + cT^ > stabilizes the EW- 
symmetric vacuum, A < causes the potential to turn over, and the 0{h^) term stabilizes 
the EW-broken vacuum. In order to obtain A < 0, we must have A < Amax^ where the 
upper bound evaluates to Amax ~ 800 GeV for rriH ~ 125 GeV. Hence, if h here is 
interpreted as exactly the Higgs direction such that niH is the Higgs mass, this class of 
models generically requires a low cutoff scale coming from trying to keep v fixed and 
A < 0. As we will see, the consequent prediction of new states at the 800 GeV scale is 
likely to be the strongest test of this class of scenarios. A potential of this form is illustrated 
in Fig. 1 . The electroweak phase transition in this effective theory was studied by [ ^-j i ]. 

As in the Class IIA scenario, the presence of the tree level barrier allows v{Tc) ~ 
V and therefore v{Tc)/Tc may be enhanced by reducing Tc. Once again using standard 

Here, we assume that the operator coefficient of the dimension six Higgs kinetic term is vanishing. 
More generally, a larger parameter space is consistent with SFOPT [49]. 



(3.31) 




(3.32) 



(3.33) 
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techniques, we calculate the phase transition temperature and the EW order parameter to 
be 




/i2 /A2A2 



Tc = \l^\l—^-l (3.34) 



v(Tc) rc~ 2 
^ ^' - I . (3.35) 



ic V -A . _ J/^ 

^ A2A2 



The optimal limits for enhancing v{Tc)/Tc are given by: 



c ^ A This limit was discussed previously in the context of the Class IIA. 

4yLt^/A^A^ — >■ 1 Using the relationships Eqs. (3.32) & (3.33), this combination of 
parameters can be expressed as 

_ 4 1 - 2AVA^,,, 
A2A2 3 (1 - AVALJ^ ■ ^ ■ ^ 

Then, the limit is obtained when A — )• Amin where Amin = ^max/ = v'^/niH- For 
tuh ~ 125 GeV this evaluates to Amin ^ 480 GeV. As we approach this limit, the 
phase transition temperature, given by Eq. (3.34), goes to zero. We found a similar 
behavior in Class IIA, and once again we can identify this degeneracy limit with an 
EdSP [10]. 

The (Why operator is able to evade the standard phenomenological constraints. 
Since it preserves the custodial SU (2) , there is no anomalous contribution to the p 
parameter, even for a low cutoff [ ]. However, if other dimension six operators are 
not forbidden, they may be constrained by electroweak precision tests. The Higgs 
cubic self-coupling, given by 

W-^(l + 2%) , (3.37) 

receives O (1) corrections in this limit. A measurement of Xhhh at the LHC is very 
difficult, but such large deviations from the SM have the potential to be measured 
with 1000 fb"^ data at 14 TeV l -]. 

A ^ We would like to take this limit A while fixing 4yU^/(A^A^) such that 
Eq. (3.35) just scales like Using the relationship Eq. (3.36), this implies 

that we must let A/Amax = const. Then, Eq. (3.32) reveals that in order to take 
A to zero we would have to take mn oc to zero. If mH is identified with the 
Higgs mass itself (recall that h in principle can be a mixture of Higgs and another 
field direction), this limit is naively at tension with the Higgs mass not being much 
smaller than the electroweak scale, but as we will see below, this is not necessarily 
a problem. 
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Figure 6. Left: The EWPT order parameter v{Tc)/Tc Eq. (3.35) (black), Higgs mass mn (blue), 
and UV cutoff A (red) as A is varied. The parameters mn and A are in units of 100 GeV. Right: The 
EWPT order parameter with niH = 125 GeV in the vicinity of the EdSP A = Amin ~ 480 GeV. 
The solid black line shows the result of a calculation using the full one-loop thermal effective 
potential whereas the dashed line shows the approximation Eq. (3.35). The Higgs cubic self- 
coupling Xhhh (green) receives 0(1) corrections in the vicinity of the EdSP. 



The non-renormalizable (WH)^ term has been studied by [50, 5 '] in the context of 
the electroweak phase transition and phenomenology. In their context, the h of Eq. (3.31) 
represents the SM Higgs field direction without mixing with another field degree of free- 
dom. We have calculated v{Tc)/Tc in the two limits (other than c 3> A) discussed above. 
First, we allow A to vary while fixing 

V/A^A^ = 0.2 (3.38) 

(a value away from the EdSP). The results, shown in Fig. 6 (left panel), indicate that 
v{Tc)/Tc grows like l/V— A as A approaches zero. For A ~ —0.15 the Higgs mass is 
consistent with the LHC signal at thh ~ 125 GeV, the phase transition is strongly first 
order, and the cutoff is low A ^ 500 GeV, which may be problematic if the LHC does not 
discover new states at this energy scale. On a positive note, as recognized by [50, 51], this 
mechanism does not rely on there being small dimensionless couplings. 

The behavior of v(Tc)/Tc nearby the EdSP is shown in Fig. 6 (right panel) where 
we have fixed uih = 125 GeV and varied A. This figure illustrates that v{Tc)/Tc grows 
as A decreases toward the EdSP where Amin ~ 480 GeV and Tc vanishes. For smaller 
values of A, electroweak symmetry breaking does not occur (and the EWSB vacuum is 
metastable). For large values of the cutoff, the Higgs self-coupling A becomes positive 
and the PT proceeds as in the SM without any enhancement. Comparing the EdSP case 
to that of Eq. (3.38), we learn that because of the EdSP enhancement of the SFOPT, the 
maximum cutoff can be taken to be somewhat larger (although limited by Ajnax)^ which 
may beneficial from a model building perspective if data pushes the possibility of BSM 
states to higher energies. The figure also shows that the Higgs cubic self-coupling Xhhh 
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receives 0(1) corrections in the neighborhood of the EdSP, however it will be difficult to 
measure this parameter at the LHC. 

3.4 Class III: Loop Driven 

In the presence of qualitatively important quantum corrections, non-polynomial field de- 
pendence may play a crucial role in rendering the electroweak phase transition strongly 
first order. For example, a competition between the terms and In h'^ may generate a 
barrier in the effective potential. Alternatively, we can say that A is positive at high scales 
and runs negative at the electroweak scale. As a prototype of this "Loop Driven" model 
class, we will consider this running quartic coupling scenario whose effective potential 
may be written as 

V^h, T) - ^ (/.^ + cT^) + h' + lhH^^^. (3.39) 
The parameters /i^ and A may be exchanged for the Higgs VEV v and mass^^ niH using 

fi'' = -'^ + Kv\ (3.41) 

The dimensionless parameter k parameterizes loop- suppressed corrections to the effective 
potential arising from interactions between the Higgs and the other fields. For example, 
in the SM one finds ksm ~ (6M^ + 3M| - l2Mf) jilQn'^v^) ^ -0.018. The loop 
induced term can help provide a barrier, as shown in Fig. 1, if /x^ > stabilizes /i = 0, 
A < turns the potential over, and k > stabilizes h = v. To allow k > 0, the BSM 
physics contributions should be dominated by bosonic fields, since fermion loops bring in 
an additional minus sign. The generation of an SFOPT through large loop corrections was 
studied in a singlet extension of the SM by [8] and in the Two Higgs Doublet Model by 
[53, 54]. 

The calculation of the EW order parameter from Eq. (3.39), requires the introduction 
of special functions (due to the non-polynomial field dependence). A more transparent set 
of expressions is obtained by performing an expansion in e = 1 — Kv'^/rnj^, which we will 
see is a small quantity in the region of parameter space that turns out to be favorable for 

'^Since the loop contributions are important in this model class, we must be careful to distinguish the 
parameter niH, defined as niH = \/V^^{v, T — 0), from the Higgs pole mass. They differ by a correction 
that depends on the renormalization conditions. Since we are primarily interested in the parametric scaling 
behavior and not numerical precision, we use uih to characterize the mass scale of fluctuations about h — v 
and implement LHC Higgs data by setting niH ~ 125 GeV. 
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SFOPT. Doing so we obtain 



ITT-H ^ A , 1 , 37 2 



„(TJ„2„V?1 ^i_?,_i03^. + . , (3.43) 



Tc rriH Ve \ 8 384 
The optimal limits for enhancing v{Tc)/Tc are given by (recall e = 1 — nv'^/mj^) 



K,v^/m^ — >■ 1 In this limit, the quantum corrections are large, i.e., k n^ax = rnj^/v"^ ~ 
0.26 for niH = 125 GeV. Since k contains a suppression factor of l/lGvr^, obtain- 
ing K = O (1) requires either many additional (bosonic) degrees of freedom and/or 
large couplings to the Higgs. This limit is then bounded by perturbativity constraints. 
Moreover, the large loops which generate k may also contribute to Higgs production 
and/or decay processes. We discuss this scenario further in Sec. 4 in the context 
of Higgs diphoton decay. Finally, we can once again identify this limit as an EdSP 
in which Tc vanishes and the EW-broken and EW-symmetric vacua are degenerate. 
Above K = Kmax clcctroweak symmetry breaking does not occur. 

"^fj <IC v^/c This limit is excluded in light of the Higgs discovery. 



As an example of a model in the Loop Driven class, we will discuss a singlet extension 
of the SM presented in [ ]. The SM Lagrangian is extended by 



N N 



AC = J2 (dSif - CH^H Sf (3.44) 

i=l 1=1 

where the real, scalar fields Si are singlets under the SM gauge group. We assume 
that > and the Si do not acquire VEVs. Instead, they modify the electroweak phase 
transition by radiatively generating a correction to the effective potential which is given by 



647r2 



In 



C^_3 

g2 2 



(3.45) 



when renormalized in the MS scheme at the scale Q. This term can be matched onto the 
logarithmic term in Eq. (3.39) by choosing k = A^CV^Stt^ and = Q^C^ exp [3/2]. 
With this identification, the limit in which nv'^/m'jj approaches unity corresponds to ( 



approaching Cmax = 2y^7rmH/vN which evaluates to Cmax ~ 2.5A^ ^^"^ for mn ^ 
125 GeV. 

Choosing mn = 125 GeV, A^ = 12, and Q = mt = 172 GeV, we calculate v{Tc)/Tc 
using the full one-loop thermal effective potential, and we present the results in Fig. 7. As 
expected, v{Tc)/Tc grows upon approaching the EdSP where kv"^ /rnjj- = 1 corresponds to 
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Higgs Portal Coupling ( ^ ) 

Figure 7. The EW order parameter evaluated 1) using the approximation Eq. (3.43) (dashed), 2) 
using the toy model potential Eq. (3.39) but without any further approximation (dotted), and 3) 
using the full one-loop thermal effective potential, as described in the text (solid). All three calcu- 
lations reveal that v(Tc)/Tc grows upon approaching the EdSP at (max ~ 1.36. For a sufficiently 
low cutoff (A ~ 1 — 10 TeV), perturbativity is maintained up to C ~ 1-5 [ ]. 

Cmax ~ 1.36. For larger values of (, electroweak symmetry breaking does not occur. For 
sufficiently small values of the PT becomes SM-like and no longer strongly first order. 
The discrepancy between the approximations and the full one-loop calculation of v{Tc)/Tc 
can be attributed to the implicit use of the high-temperature expansion in Eq. (3.39) and 
setting c = csm ~ 0.36 without accounting for the S field. The contribution from S 
is suppressed at large ( (where the approximation agrees well), because S is heavy and 
its thermal contribution is Boltzmann suppressed. At smaller values of the S field 
effectively renders c > csm which tends to increase v{Tc)/Tc, as indicated by Eq. (3.43), 
and confirmed by Fig. 7. In this model, the additional singlet scalars will not have an 
appreciable impact on collider physics. We discuss the more general case of charged or 
colored scalars in Sec. 4. 



4 Diphoton Excess and SFOPT in the Higgs Portal 

One tantalizing hint of new physics in the recent LHC announcement is the observed 
excess of events in the final states with two photons. The 77 final state, which is associated 
with Higgs production by gluon fusion, is observed at a rate that exceeds the SM prediction 
by a factor of approximately 1.5, while the 77jj final state is enhanced by a factor of 
approximately three [24]. Although not statistically significant yet, fits to the entire data 
set seem to favor an enhancement of the diphoton decay rate T{h — )■ 77) by a factor of 
approximately 2 — 3 with respect to the SM prediction, as well as a suppression of the 
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gluon fusion production cross section cr{gg — )■ /i) by a factor of approximately 0.5 — 0.6 
(see, e.g., [24, J9] and references therein). 

Since gg ^ h and /i — )• 77 are both loop-induced processes in the SM, these channels 
are particularly sensitive to new physics. For instance, the appropriate enhancement and 
suppression can be achieved by letting the Higgs couple to a new scalar S via the Higgs 
portal [42, 55, 56]. If S is charged, graphs containing an S loop will contribute to the 
amplitude for /i — )■ 77 and interfere with the t— and VT— loops that dominant the SM 
contribution. Generally, a negative value of the Higgs portal coupling is favored if the 
/i — > 77 rate is enhanced, because then the S-loop will interfere constructively with the 
SM contribution. Furthermore, if S is colored it will also interfere destructively with 
the SM gg ^ h. As we have seen, the Higgs portal operator also provides a means of 
rendering the electroweak phase transition strongly first order. It is then interesting to ask 
whether the region of parameter space that can accommodate a SFOPT can also allow for 
enhanced diphoton decay in such simple models where a single Higgs portal operator is 
responsible for both phenomena. We will see that generically, the SFOPT condition favors 
a positive value of the Higgs portal coupling and, therefore, is at tension with the diphoton 
enhancement in such minimal model settings. 

In order to demonstrate that SFOPT favors a positive Higgs portal coupling, let us 
consider such an interaction between the Higgs and a scalar field S, as given by the La- 
grangian 

-CD filS*S + 2XH^HS*S . (4.1) 

The phase transition calculation is independent of the quantum numbers of S at the one- 
loop order, but instead only depends upon the coupling of S to the Higgs. However, in 
order to obtain an enhanced diphoton decay rate, we need S to carry an electric charge. 
Consequently, we must ensure that S does not acquire a VEV.^^ In that case, the field- 
dependent squared mass of the 5* field is given by 

mls,s{h, T) = /i| + Xh^ + Us{T) , (4.2) 

where 115 (T) is the thermal self-energy correction. In the appropriate limits, this simple 
extension of the SM Lagrangian can yield any one of the phase transition models classes 
discussed above. These are as follows: 

Class I. Thermally (BEC) Driven. The BEC term receives a contribution + 
Xh"^ + UsiT))^^"^. As discussed in Sec. 3.1, we must tune /i| ^ —Ils{Tc). However, 

'^If S is colored, then the two-loop contribution from gluons can have an appreciable impact on the order 
of the phase transition. This is, for example, the case in the MSSM 

'""This discussion presumes that S* is a singlet under weak isospin. More generally, the electrically neutral 
component of S may acquire a VEV without breaking U(l)^^^. However, unless this VEV is much less than 
V, it will be at tension with electroweak precision measurements. 
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in this limit the mass of the S field is m^j^ g(v, 0) = —Ils{Tc) + Xv'^. We cannot let 
A < 0, because this would render S tachyonic and induce a VEV. 

Class IIA. Tree-Level (Ren. Ops.) Driven. Since S cannot acquire a VEV, the only 
way in which tree-level terms can enhance the strength of the phase transition are 
if S had a nonzero expectation value in the early universe which returned to zero 
during the electroweak phase transition. This scenario is realized by letting |u| < 
such that S obtains a nonzero expectation value in the early universe, but ensuring 
— Af ^ < ;u| < such that S has a vanishing VEV today. Once again, we find that 
A > is required for a SFOPT in this model class as well. 

Class IIB. Tree-Level (Non-Ren. Ops.) Driven. The non-renormalizable operator 
Hy / A"^ may be generated by integrating out the field S. The leading order 
contribution to this operator coefficient is proportional to (l/167r^)(A'^/Af|). Since 
this model class relies upon (WH)^ having a positive coefficient in order to stabihze 
the potential against a runaway direction, we must take A > 0. 

Class III. Loop Driven. This model class relies upon the addition of a term to the 
effective potential that goes like In h'^ and its competition with the h'^ term to 
generate a barrier in the effective potential. The Higgs portal operators Eq. (4.1) 
will instead generate a term of the form h^ln(fXg + A/i^). Unless <^ |Af ^|, this 
term will simply scale like h'^ and there will be no competition between terms and 
no barrier. However, if A < 0, then in this limit the S field develops a tachyonic 
instability and acquires a VEV. 

This analysis may seem to suggest that A > is generally favored by SFOPT. How- 
ever, this is not the case. If we were not interested in enhancing Higgs diphoton decay, 
then we could achieve a SFOPT by coupling the Higgs to a singlet scalar field using the 
operators Eq. (4.1). Choosing A < 0, there exist models in the Tree-Level (Ren. Ops.) 
Driven class which achieve an SFOPT when the singlet has a VEV, which is restricted by 
Higgs mixing constriants (see, e.g., [ ]). Furthermore, in a non-minimal model setting in 
which we introduce additional singlet and charged scalars, the singlet(s) can enhance the 
phase transition while the charged field(s) can enhance the diphoton rates. 

5 Conclusion 

In this paper we have proposed a classification of the electroweak symmetry breaking sec- 
tor which may yield a strongly first order phase transition - a necessary ingredient for 
electroweak baryogenesis. For each model class, we assumed that the last phase transition 
associated with the electroweak symmetry breaking sector was an electroweak symmetry 
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breaking transition (i.e., no broken vacuum to broken vacuum transition), and we inves- 
tigated the impact of the data that is currently available from the LHC: 1) the discovery 
of a 125 GeV Higgs-like scalar, 2) the absence of a large exotic (e.g., invisible) decay 
width, and 3) the absence of a universal suppression, which would indicate mixing be- 
tween the Higgs and a hidden sector scalar field. We find that the mass measurement 
severely constraints models (such as the MSSM [14, 15]) which drive a strongly first order 
phase transition with thermal loop effects. The invisible decay and mixing constraints are 
at tension with models which rely on light singlets coupled to the Higgs. 

One recurring theme of our analysis is the ubiquity of enhanced symmetry points. We 
find that the "optimal" limit for SFOPT often corresponds to a parameter point at which 
the symmetry group of the theory is extended. In the case that the group is enlarged by 
a continuous symmetry, either the Higgs mass constraint or the exotic decay and mixing 
constraints will come into play. The case of a discrete symmetry is less restricted [10]. 

We have also discussed the possibility of employing the same Higgs portal operator 
to both render the EWPT strongly first order and to account for the diphoton excess ob- 
served by ATLAS and CMS. We find that these two goals are at odds with one another 
in the minimal model setting: the phase transition favors a positive Higgs portal coupling 
whereas the diphoton enhancement favors a negative coupling. A model which can ac- 
commodate EWBG as well as fit the LHC data will most likely require two distinct new 
physics operators. However, it is worth noting that the diphoton excess does not have a 
great statistical significance, and the data remains consistent with the SM at the 75 % CL 
[22] or approximately 2a [23, 24]. It is still entirely possible that the particle recently 
discovered by ATLAS and CMS is the SM Higgs [j / ]. 
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